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SOQ USER GUIDE UPDATES 


June 1980 Updates to SO0Q80128 


INTRODUCTION 


This document defines the changes made to the SOQ code (S0Q80128) 
between January and June of 1980. The changes either correct short- 
comings found in the code or, more usually, document the increased 
capability being continually built into the code. The SOQ code i 
raintained as S0Q80128 June PL,ID = ARLOJRA as a NOS/BE-1 CDC update 
format file. 


UPDATES 

1. *ID FIXZRN 

This update redefines the coeffiencets to be input to the Zernike 
subroutine. This new convention is more physically meaningful in that, 
at least for lower orders, the coefficients are in waves. For example, 


to impose one wave peak to peak of defocus (P,) on a beam, one would 
input P(4)=1. The phase applied is now: 


o(1,J) = pp mz (14d) 


The subroutine affected is ZERN. This update does not effect the rest 
of the code. 


2. ¥*ID FIXJTR 

This update ensures a correct definition of DF in subroutine JITRBG 
since when JITRBG is called from subroutine QUAL, the X-coordinate 
array contains RA/D coordinates, not the spatial coordinates. 

Only one line of the code is affected by this update. 


3. *ID ROTZRN 


Due to different coordinate system orientations for data, it became 
necessary to allow for this variation within subroutine ZERN. 


Define the data x and y coordinates to be XROT and YROT, and the SOQ 
x and y coordinates to be XIN and YIN. The rotation angle is then 
defined to be 6 (in radians). 
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COSROT = COS(@) 
SINROT = SIN(0) 


XROT XIN x COSROT 


+ YIN x SINROT 


YROT ~XIN x SINROT + YIN x COSROT 


Application of Zernike polynomials to and SOQ point located at 


(XIN, YIN) would then be 
bility of axis flips are 
FLIPX or FLIPY not equal 
include FLIPX, FLIPY and 


calculated using Z(XROT, YROT). The possi- 
also accounted for and are flagged by 

to zero. Namelist ZERNS is modified to 

the rotation angle (in degrees) ZTHETA. No 


common was modified. This update modified only subroutines GDL and 


ZERN. 
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SECTION V APPENDICES 
APPENDIX A 
SOQ USERS GUIDE UPDATES 
JANUARY 1977 TO JUNE 1978 


INTRODUCTION 


This appendix documents those changes made to the initial SOQ code between 
January 1977 and June 1978. The changes incorporated in the code are those that 
have become generally useful for the physical optics simulation problems which 
have been solved using the SOQ code. The Users Guide Updates are also prepared 
to clarify and correct the initial description of the SOQ code, as documented 
and delivered to AFWL on 1 March 1978, in the Preliminary SOQ Users Guide. 

This document supersedes previous written material on SOQ code documentation. 


The organization of the SOQ Users Guide Updates is 


SECTION AI New Subroutines 
1. Theory 
2. FORTRAN Updates 
SECTION AIT Code Changes/Corrections 
1. Theory/Reason for Correction 
2. FORTRAN Updates 
SECTION AIII Users Guide Corrections 


SECTION AIV SOQ Code Access 


AI. NEW SUBROUTINES 
1. THEORY 


a. Beam jitter -- Relative motion between optical elements, such as 
mirrors, induces time varying positional displacement of the optical field. 
The typical term for this phenomenon is beam jitter, and the principle effect 
is to broaden the sime-averaged effective beam illumination area, while 


reducing the time averaged intensity. 


Beam jitter is both a near field and far field concern. Jitter in optical 


trains can overload apertures or cause energy deposition on areas outside the 


normal propagation path as well as cause a deterioration in the peak on-axis 


irradiance and integrated spot power. 
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The time-averaged effect of beam jitter may be modeled as the convolution 


of the intensity profile with an appropriate probability density function (PDF) 


for the jitter statistics. The current SOQ model assumes that the jitter PDF 


is Gaussian with known mean and variance. The model allows the user to speciry 


the Gaussian parameters and, for the selected beam jitter analysis location, to 


determine the near field and/or far field effect of beam jitter. 


The following is a brief description of the analytical and SOQ Fortran 


implementation of beam jitter calculations: 


b, Relevant formalism -- The effect on the beam may be found by convolu- 


tion of the Gaussian jitter probability density function with the SOQ predicted 


intensity distribution: 


J 
1’ (x,y) - ff: (hy) J (x-x’, y-y’ ) dx’ dy’ (Al) 
-00 


The 1-D Fourier transform of the Gaussian function looks like: 


a peer 


H(f) 


0 


FORTRAN UPDATES 


Each assumes that the jitter 


The jitter model can be called in two ways. 
variance is the product of a jitter angle and the propagation distance from the 


jitter source. 
(A2) 


= Jitter angle (lo, in microradians) 


= Distance from jitter source (in cm) 


When the far field model is called from QUAL, the jitter angle has been 
incorporated into namelist QLOT while the propagation distance is the focal 
length found in QUAL. The jittered intensity is returned to array CU as a 
phaseless field so it can be plotted, or written to a permanent file. 


The other method of activating the jitter model is to call the near field 
jitter model from GDL with IFLOW = 23. For this model both angle and jitter 
distance are entered in namelist JITTER. 

Namelists modified: 
Far Field 
QLOT: SIGANG (rad) is added to specify the jitter angle 


Near Field 
JITTER: Contains - 
JITANG (urad) 
JITDIS - Jitter distance 


ATI. CODE CHANGES/CORRECTIONS 


THEORY/REASON FOR CORRECTIONS 


a. Bare resonator calculations -- The SOQ resonator/optical-train cal- 
culation code may be used to simulate, in Cartesian coordinates, bare reso- 
nators. This added option is frequently used in the initial simulation studies 
of a resonator or a class of resonators. 


The bare resonator optical configuration may be compared to its geometric 


counterpart using the SOQ code by simply invoking the IBARE option and asso- 


ciated updates now contained in the fundamental code. The fundamental approach 


in bare resonator calculations on the SOQ code is to allow the user to use the 
same input and code for bare, semibare and loaded cavity calculations. Various 
options under the bare cavity calculations have been incorporated and are now 
described as input values for IBARE in Namelist START. 
IBARE = 0 (Default) 
Loaded cavity calculations are performed as usual following the stan- 
dard input which the user has supplied. 
IBARE = 1 
Using the same input, the user will now perform bare cavity calcula- 
tions in which the resonator is normalized to 1W of circulating power. 
Mirrors are defaulted to have 100 percent reflectivity, and no power 
dependent or flux dependent distortion. The SOQ output is modified 
“to printout the resonator eigenvalue. 


IBARE = 2 
Semibare resonator calculations are performed in which the user can 
perform bare resonator calculations that include optical aberrations 
generated by a flowing saturable gain medium. These aberrations may 
strongly effect mode shape/phase. This option provides a convenient 
method of studying their perturbational effect on the bare cavity mode. 


For the semibare option, an additional update has been included in which 
the namelist MIRROR user may specify the desired power at each of the resonator 
mirrors. This allows the user of the semibare updates to apply mirror distor- 
tions as though the bare cavity mode had a significant power level. Specifica- 
tion of the DESIPW value in namelist MIRROR to some value other than 0.0 will 
cause the field incident on the mirror to be scaled to that power specified by 


the numerical value of DESIPW. Appropriate mirror distortions will be applied 
at the desired power level. The field leaving the mirror will then be rescaled 
to its incident power level. Subsequent calculations are done as specificed by 


the user typical namelist input. 


An additional variation is allowed in which the parameter FLAG of cavity 
input namelist CAVTY2 can be used to execute a resonator with loaded gain, but 
no fixed phase perturbation in cavity. The input would correspond to FLAG = 0; 
IBARE = 0. Usual loaded resonator calculations are performed with mirror dis- 
tortions as specified by the user. 


All of the above variations of cavity/resonator calculations may be run 


from the standard loaded cavity input. 


At the rear of this section are Fortran listings of the code updates which 
have been included in the basic Cycle III SOQ code previously documented. 


b. Mirror non-normal incidence angle -- In many optical train calcula- 
tions the propagating optical field is incident on the mirrors in a nonnormal 
manner. Since, in general, the mirror surface may have a spherical figure, the 


field leaving the mirror will exhibit phase front aberrations introduced by non- 


normal incidence of the field on the curved surface. 


The SOQ MIRROR subroutine has been modified to incorporate an astigmatic 
aberration due to the nonnormal incidence on a spherical surface. The follow- 
ing is a brief description of the generation of the astigmatic aberration 
applied. 


Astigmatism in Resonator: 


General astigmatism is introduced when a wavefront is incident on a spheri- 
cal (parabolic) surface in a nonnormal manner. This aberration occurs at each 
spherically-distorted turning flat, for example. 


z 
y x - Sagittal Plane 
y - Tangential Plane ‘O) 
dg Bl.) GROSS 
Ss 37, : Ro 4 = Incident angle 
R, S #2 Object distance 
imi Cs Qe 
es 2 cos¢ R_ = Mirror surface curvature 
. Cc 
(spherical) 
Ae 2 ot x? S, = Sagittal plane effective 
SOQ, r 285, curvature 


Thus 28° is the resultant phase curvature being imposed on the beam. A 
cylindrical mirror can be used to model this with 
F R. - 2s° = -V2R, (neg since Ro is convex) for 9 = 45 
: (A3) 
Therefore, to represent the astigmatism introduced in the x-plane by a 
spherically-distorted turning flat, a cylindrical mirror is employed with a 


radius of curvature 


a -v2 R, 
Ss 
Re is the power induced radius of 
curvature which is input or deter- 
: mined by the SOQ code. 
« Similarly, the tangential nlane is described by 


46 . an (_y? 
SOQ, iN 2s’, 
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-Rcos¢ 
ago (A4) 
= —& for o= 45° 
2v7 
The new mirror subroutine including astigmatic effects has the form 
2r x? 
8509 7 | 3er * Ser 
s Tr, 
R R_ cos 
Ss’ = —_o__ Si. = c ? 
s 2 cos ¢ rT 2 (AS) 


The only additional input change is to the MIRROR routine namelist which 
is expanded to include the variable PHIAST, the beam incidence angle in degrees 
(default is PHIAST 2 0). 


c. Beam rotation -- The mirror model has been updated to describe beam 
rotation introduced by optical elements. which are oriented in a skewed fashion. 
Many examples of this type of orientation are encountered in resonators and 
optical trains. The principle effect of skewed, or out-of-plane, orientation 
is to convolve or smooth the mirror distortion-induced aberrations over the 
total number of optical elements. 


Rotation of the beam is accomplished by analytically rotating the mirror 
with respect to the beam, rather than rotating the beam within the mesh and 
then applying the mirror. By rotating the mirror with respect to the beam 
two modeling advantages result: First, analytical rotation of the mirror with 
respect to the beam is accomplished with no interpolation loss of information. 
Second, since the rotation is analytical, computer time is saved by not having 


to evaluate the field numerically. The following describes the rotation 


equations used in the code. The following sketch shows a base and rotated system. 


y 


a = Rotation Angle 


xX = X cosa + y sina (A6) 


y = -x sina + y cosa (A7) 
Then, 
2 
2. 
q Here, 


(x,y) are the SOQ coordinates 


(x,y) are the transformed (rotated) coordinates 


The SOQ field is modified as 


‘ A9 
cu A exp (49) CULy (a9) 


OUT 


where A represents the completed transmittance effects included in mirror. 


The variable added to the SOQ MIRROR namelist input is PHIROT, which is the 
beam rotation angle in degrees. The default value is PHIROT = 0.0. 


2. FORTRAN UPDATES 


The attached printouts contain a listing of the updates, denoted as ROT, 
used to effect these changes. 


Litheal wt 


@LOENT wat 
“DELETE ClOaste.) 
ATOP (Dee) oXSCA 6G) ABC (1402909) o FI TLEI (20) eXUMANN (4) 6 
*ORLETE CLOFL A. 
DIMPNSTUN TN [9(5 626) o Tail (99) oAHC (14020 09) eCFFL (16386) e LUSK (409) 
®URLE TE CITOASTGL2eCILOASTG.3 
DATA WANES eDOUTYSUINY oPH(AST ePHIHUr /500./ 
PUVELETE C1LOASTR.G 
A DELTA SISTR s2OUTY SULNY eH ANULS ePHIAST PMIROT NE SI PW 
®INSERT C1NaSTSG.5 
Cc MHTMOT TS Tee HEAM ROTATLI Ul ANOLE AT TraT STATI One DEG 
¢C NESTPwW TS TmE NESTRED POwEN LEVEL AT TeaT STATIUN 
*®INGEWT C1NASTG.A 
AAC (1 %eIMTRee) ON[POT 
. AAC (Le@eIMIRe?) 2 NEST Pw 
2e OESTROw 2 AHC(LEes{IMIRoeg} 
HELE TE ClOasSTGed 
; X OTSTE eWANULS ePYOUT ORY INGMHTAST ePMINT NEST Pw) 
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PINSERT MIRWOR, 2A 
PHTROTsPrInes 
*PINSERT CiNasTsu.t2 
PHIRURS@PHIMUTOEMIT/ lL BU.e 
PHIROTZ0.0 
aR ITE (6086) PHI RU 
SINPRASIN(PRTROR) 
COSPRSsCUS (PMT RUR) 
*NELETE C10aSfG.15 
XPRMBX (J) OCOSPHOXK (1) OSINPR 
YPRMsek (J) OSINPROK (1) @CUSPH 
PHASE SAK Y®( (XOERM2O2/HMSAG) © ( YRRMOO2/OMTAN) ) AK YONELL 
*VELETF C1HASTG.26 
APRMax (J) OCOSPHOX (1) OSINPR 
YPRMaek (4) OSINPWOk (1) PCUSPH 
PHASE SAKY® ( (KPRM®O2/HMSAG) © ( YPHM@OCs/OMTAN) } 
®DELETE MIRROR .BS 
PHI Re (PHIASTOEPT)/180. 
RMSAGSRUC/COS(PHIR) 
RMT ANSHOC COS (PHIR) 
PHI ROR SePHIROTORI/1 40. 
PHlhOTsO.0 
SINPRSSIN(PHIRUR} 
COSPRSCOS (PHI RUN) 
*PNELETE MIPROR.D1 
XPRMEK (J) OCOSPREK(T) OSI NER 
YPRMSOk (J) OSINPROK (1) @COSPH 
PHIMIR ZAR Y® ( (XPRM@O2/RMSAG) © ( YPWM®O2/QMTAN) } 
CINSERT MIPROR.108 
WHITE (6046) PHIROR 
WRITE (46421) RMSAGeRMTAN 
-J-) FORMATI(2040@ MIRROR RUTATIUN = %9612.50°¢R 00S) 


ATIT. USER'S GUIDE CORRECTIONS 
SUBROUTINE FUHS 


a. Purpose -- Subroutine FUHS is used to calculate the phase change 
due to heat release as the molecules in the lower laser level decay to the 
ground state. The FUHS modeling includes the assumption generally made for 
supersonic flow and assumes the heat release has only a small perturbative 
effect on the flow. 


b. Formulation -- The equations used here are based on those described 
by Biblarz and Fuhs, (Ref. 10) and-Fuhs, (Ref. 11). 


The usual continuity, momentum, and energy equations for steady flow 
with heat addition are used as the basis for the analysis: 


Continuity: V- (pu) = 0 


Do = 
Energy: Vepu h + & 


Momentum: 


= 


These are linearized assuming 


* ) , A» 
G2ei(+1u)+ijvy 


Resulting in 


Continuity: 9, uw, + vy + Up’ = 0 (ALL) 


rey 
8 
rs} 


(Ys = uw; etc) 


Momentum: 


Energy: 


The solution to these equations was first shown by Tsien and Bieloch, 
(Ref. 12), resulting in the following equations for a heat source q in super-~ 


sonic heat addition. 


re -(y - 1)q 46(x - 8y) 
2yp8 


(A15) 


». &-}) 6(x - By) (A16) 
iy cae 
p= (r= 1)Mg 8(x - By) | (y= a Sly) I(x) (A17) 


2a°B a“U 


Defines a Mach line 


Speed of Sound 


For a small volume, the heat addition is q = h{x,y) dx dy. The effects of all 


sources are then added; for example, 


Pies Gas) é 
pt ffrceor 5(x-gy) dx dy 


s 


(7-1) ; 
=- “S54 f h(x=$y) sin yds (A18) 
rf 


where the integral is taken along a streamline (x = y) and siny = 1/M. 


s is related to x and y by 


x = S$ cos# y = s sing (A19) 
By the Faltung theorem, for Fourier transforms, this can be written 
I’ (x,y)= P| F [t¢x,¥) | °F [sen] (A20) 


The Fourier transform of the intensity is performed by the FFT, while the 


Fourier Transform for Gaussian density functions can be found analytically 


as 
229 
" a exp poy exp [ ~2n0" (£7 + £ | (A21) 


From sampling theory, the discrete values of f. and fy can be found since 


A£ = + (A22) 


where 


L is the width of the SOQ calculation region (DCALC) 


£,() is then (I-1)4f. Recall from the discussion in FOURT, the DC value is 
returned in the first position and the last half of the transformed data are 
really negative frequency information shifted by one period, illustrated 


below in one dimension. 


Ax = Sampling rate in real space 


The equation for density change is, therefore, 
s 


ik E : os { h(x,y) lan? sin u ds 


? 2a°8 


(rs off dy’h(xyy’) 8(y = y!) T(x - x) 
aU 


The first term describes the compression waves along the streamlines due 
to heat addition, while the second describes the wake resulting from those 
compression waves. 


The heat release h(x,y) for a laser can be written: 


% ~ xX) 
h(x,y) = c feces Rt 


NEP 


where t is the time constant for the depopulation of the lower laser level. 
If the depopulation were instantaneous (t + 0), the heat release would be 
proportional to the intensity, since every molecule emitting a photon would 
then immediately relax to the ground state with an accompanying increase in 
translational energy. It has been shown that the above equation for the 
heat release can be used in all regions of the cavity with only small error 
introduced. 


: The constant c can be found by conservation of energy as shown 


following. 


Consider the following 3-level molecule: 


E2 


Ey 


Eo 


The quantum efficiency n is defined as the ratio of the energy out divided 
by the total energy available, so for the gain/phase segment under 


consideration. 


(no. of molecules) (E, - E,) AP (A25) 


ks (no. of molecules) (E, - Ey) * aH + ap 


| 

4 
assume, for this calculation, that (0,0) is at the corner of the sidewall 
! fi 


Where 
AH = (no. of molecules) (Ey - E,) ° 
the above expression can be inverted to give 
AH -(=} (A26) 
with 
AP =fforcey dx’ dy’ 
and 
wave frcar) dx’ dy’ 
and the NEP. Then 
PP . -(x - x’)/Urt 
OH = cAz oy fa [ sncey Je dx’ 
e e e 


-(x - x’)/Ut 
T(x - x’ )Al(xjy)e ) dx’ (A27) 


~e 
> 
——, 
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Where, recall 


| x > x! 


Since the numerical kinetics routine returns information about the wake 


I(x - x! 
“Yo, x < x’ 
So, 
- FP P “(x - x!)/Ut 
AH = cAz fe eke I(x - x’Je 
rp - -x” /Ut 
= cAz f dy [« AL (x3 : dx"e 
7 «? 
= chz J» fo AI (x4 y)} “(che (A28) 
Or, 
AH = c(Az) UTAP 
& 
re oe; 
« (lca n 1 i 
| . ( 7 ] Urdz t 
{ 


region itself and not just the heat addition terms, this information must 
be the data used. Thus, for the analytical kinetics model, one must find 


woe ce 


the value for the wake integral: 


x 


4 - -(x!l- x") /Ut 
w(x.y) = fax'ncty = c fexfexrarcer.yre 
eo e eo 


2 cf 20 - x) factee- x") AI(x”,y)e 


-(x! = x") /Ur 


~ ~ =(x! = x") /Ut 


= z femrarcery fax’ tex - xf) I(xy’- x%e 


~ -(x = x") /Ut 
2 c fextarce» I(x - x”) Jdx’e (A29) 


x 


Stina ns Sie ae a i fe 


So, 
x 
-(x - x") /Ur ‘ 
w(x,y) = c | dx"Al(x’,y)UT(1 - e ) (A30) 
° 
So, recalling 
& epkeeomy se (A31) 
n UtAz 
And 
1-G from SIMPGG 
ia -(x-x") /Ut 
2 {1 -G\f/l-n PPD (x4y)(l-e ) 
wey = (Le $)(Len) [ 
Az \l G n : (A32) 
Now both numerical and analytical kinetics models return the same array, ; 


namely the value of the wake integral throughout the cavity. The effect of 
heat release due to lower level depopulation can be calculated without 
regard to the particular kinetics model chosen. The Fuhs effect is calculated 


in the following manner: 


x (T) 
1) }-w(x(I-1)) 
H(I,J) = a freee = “se a (A33) 
X(I-1) 


Given this average heat release function, the integral along a 
characteristic can be performed. Note that reflection off the sidewalls 


must be included as can be seen in the following diagram: 


1 » NEP 
N 
| ‘ 
\ | \ 
\ | \ 
| /\" 
4 
¢ 
ec 
The contribution at iat due to reflection is therefore found by finding 
the total heat released along the characteristic that reflects at P,, 
then adding this to that found along PoPy: 
! | 
The phase shift is found using the Gladstone-Dale law 
n2zil +*Co (A354) i 
The phase change d¢ is : 
2" (CA 
2m (An) (dz — (=e AS * 
do = oe fan) (42) = 5 ic 0% (A35) 


This is then added to that of the unloaded density field to establish the total 
phase change at the gain/phase segment. 
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AIV. SOQ CODE ACCESS 


1. SOURCE CODE 


The following listing represents the source code necessary to update 
the SOQ to include the corrections and modifications described on the 


preceding pages. 


* Toe eT Mowh 


STOENT sake 
YINSERT GUL.S3] 
ee Tr CIRAQE .NELL) GO TO 850 
‘ WMIR s 1.0 
OELTA 0.” 
NTSTF = 0, 
AM ITTE (Re ABN) [RARE 
BOO FORMAT /2X e LOHPFO2SOHSEO.4 THAne S%e[ Deo NELTA ANDSe 
C et HISte Ser TA ZFWO AM) wMT 2 SET TU Le Pel Oneseeceaseuys) { 
ssu CONTINUE 
PINSEWT APO OK,?2D 
wtITE(7) (CUCLTZ7) eT 23] 6NUn) 


ay 


REwING 7? 
DELETE ANL e7060GOL.795 
PIELER TE BOL eA9N eudL eT 


TEC TCAV EV.) [HAKE SL 
IF (LRABE.FQ.0? 60 TQ AY] 
WRITE (29) (CFFL CLL) e TL) sind) ok 
REWIND 29 
NO GAY T21eN03 
oad Cucl) = CFFLCT) 


Mi) GAGT LaleNPTS % 
4adl K(f) = KKCT) 

OTAMH = ARCllelel) 2-5 

CALL APRTM (DTA He VesVe ede ee ede) 

POwsn, 

REAN (29) (CFFL (TZ) eo L221 eNUH) 0X 

REWIND 29 

110 68m [el eNOs 
6gR POW = Pilw 6© Cult) * COwGI(CUII)) 


Pow = POW 9(XK(2) | AN(L)) 292° ( NPI S/NBY) 
WRITE (Ae6H7) POW 
oa? FORMAT(SXy® come PUWER IN FEEDRACK NORMED TO UNITY HY Ge 
t) KE LS e7e/) 
SQTPOW = SQRT(KOw) . 
OU 684 [31 eNO 
5d CFFL (I) = CFFLIIT) / SUTRPOW 
: ®LELETE GIL A22eb0L 423 
TF (T8 ARE LEN.0)} GO TO LON? 
TF CTT TO AMIIY, NOT.RESTRII GU TO A? 
Cc ®2OV8o CALCULATE FILENVALIE PISSGHOOBZS 
TF (POwKkK .GT.9200L) PRPwR = 40} 
FIG = SART( 1.210009. She wWK) 
wr ITE (HeROIF IG 
86 FORMAT (2NKe® ESTIMATED EIGENVALUE = ®eE1S.72 0/7) 
G&) TO 3093 | 
of WRITE (6 6AR) | 
ay] FORMAT (20xe9F [AST PASSeIweUT PUwEM NUT UNITY eEIG NUT ESTD4%, 7) ; 
G) TO 1902 | 
| 
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1V02 CUNTINVE 
CALL REGAININCT oN TER) 
1003 TFCICEK.EN.0) GU TU 964 
FUeLETE GIL R828 «GOL e434 
PORLETE CAVITY.3 
wR 2LIeZL0el Aare) 


PINGERT CAVITY.141 


IF (LRARE NE .O} Gu TO Las 


PUELETE CAVITY.1Aa3 


199 WRITE (C7) (CUCT™’) ofLZsleMuT) 


PORLETE CORR] .a? 


CALL DENSY (FLAG eARU oALEN © YLENSUCZ NAA ONY AG] o [No WNSYMe THARE) 


PINSERT SON77CY1.19 


IF (1IRORE.GT.O) GO TO Le 


PINSERT CAVITY.217 
6G) TO 11 


12 CSCL7) sCMPLA(COS (PHIM) eSIN(PHI™) ) 


PUELETE NENSY.2 


SUBROUTINE DENSY (FLAG eRHU eXLENe YLENGZSt AR ONPX oNP Yo LF oI NeNS YM, 


xX TAARE) 
@LVSERT DENSY.105 
TE (CTRAWF LEQ. 1) RETURN 
IF (L46.G67T.0) GO TO le 
wRITE (60) 3) 


17 FORMAT S/1 9X eSHOee = 6 PFLAG 


4a 19X- 
RETURN 
12 CONTINUE 
PINSEST LROPI.2 
OATA [AARE/SO/ 
PINSERT GOL .467 
Nara NESTRPWsN./ 
*INSEXT LROPI.L 
CUMMON /HARES/ [RARE 
PINSERT GOL.22 
COMMON /AARFS/ [Rake 
SNELETE LROK1L.3 
K s[SARE ePLOTS 


SFNAMAN 


INSERT GOL 2430 

TNIR(GeI CAV) & IRARE 
®INSERT GOL .437 

RolOTA (Sel Cav) 
DELETE GOL.439 


De [N VEMSY%eoH 8 a40/ 


3P0S SET TO 9.08//) 


TAARF TS FLAG FUR LUALEI)« HAREs UR SEMI@RARE CAVITY 
2 0 FOR LOADED KESUNATOR (UOEF AULT VALUE) 
= 1 FOR RARE RESUNATOR 
32 2? FOR SFMI@RAKE RESUNGTUH (UINITY GAINSDENSITY PHOSE CHANGF) 


(UNTTY GATNed PHASE CHANGF) 


X = PROPAGATING PARAMETER o[2/% THARE® *o 13) 


*NELETE GOL .446 


X RESTHTOINIR (Ge (CAV) ZL I (ITCAV) sZLOCICAV) s IOI R(SeICAV)) 


PORLE TF GOL.S30 


TF(eNUTSINIT) GU TO 22 


NESTPw = 9,0 
PINSEAYT GDL .2465 

NOK 3 NPTSenPyY 
®DELETE “TQWOR, 26 


70 TF(DESTOWLEQ L060) GY TO 360 


NOd = NPYONPTS 
NOK2 3 NURe2 
C eee FIND INCIDENT PUWER 
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Shins > 32a ene i neem tel Sri 


ew 


TORT 


Lord 


AARE 


PrwIN = 0.9 
90 355 [221 eNOR2e2 

355 PPWIN = PPwitt ¢ CUR(IZ)®¢2 «© CuR(I1Z°1) 402 
PPwK = PPwIN®(X(2) 9X61) )%*29(WRTS/NPY) /1 000. 
IF (NRE GR.EQ. 1 eOR NREG EQ ee) PP wk SPP WK /WNOWOe? 
TRANS = SQHT(OFS[PWs/PRwH) 

C #@@* SCALE THE HEAM TO THE O&SIREW) POWER, 

OO 356 [2 2 1 seNNd2 


356 CUR(TZ) = CUR(IZ) @TRANS 
WRITE (6¢4010) DEST PwePrwk o TRANS 
4010 FORMAT (/5Xe4AHTHE FIELD HAS BEEN SCALED TO DESIRED POWERS 
A AXel 2HDEST PW B9G1204/8X 0] CMPOwK 20Gl2e4/ 
X  BXel2HTRANS BeG1l204/) 


DESI Pw = TRANS 
369 TF (CARS CANX) LE cO-090LUU cAND ec ANS (ANY) ot F20.0N9010U) GO TN 7] 

*DELETE Cl10astG.?7 

3 ABC (11 e IMIR o%) eARC(L2eIMI 462) eABC(L3IeIMIR«2) NESIPwW) 

TF (DESTPWe.Fe0.0) GO 10 ae 
C @90008 SCALE THE FIELD BACK DOWN. 

9U 368 172] eNUG 
356 CU(I7) 2 CUCIZ) /DESIPW 

ART TE (604000) 

augn FORMAT (/SXe30HTHE FIELO RAS HEEN SCAI FD NOWNS) 
24 CONTINUE 


PTDENT JITTER 
SDELETE GOL.20 
K [CUT eML To fOKeIUTHe 
X ICNT2GeICNT25 eI CNT2Z6% 
KX [TMe ICE K oNCT 
YNELETE GOL.26 
XNSMM (21) eRUV IAN) ePHIA( ZO? eRCURVE (4) 0OSF (4) eo TLT (4) eI CAVZ(2y) 
*0elLeTE GNL.260 
90 173 IZEROZ1Le20 
*0QELETE GOL.314 
OO 177 L7FROz1l617 
*INSERT SO077CY1.2165 


Cc 2 23) JITTERS THE de A4 AN ANGLE ANGJIT 

Cc = 26 DUMMY = LINE 260 [5 TEMPURARILY STORED IN JITTE® TFLOW, 
Cc = 25 UUMMY © LINE 250 [TS TEMPURARILY STORED IN JITTER IFLOW, 
c = 26 OUAMY © LINE 260 [5S TEMPURARILY STORED IN JITTER TFLOW, 


®INSERT CORRA,/ 
DATA LUTReJITANGeJITUIS 7000 000020/ 
SUELETE GOL .799eS0077CY1.107 
¢ 4 les 17/7 LAL 197 20/4 el/ 22/7 234 26/ 28/ 26s 
X 015061704180 0190 0200 9210 2 36902350 02600250 6260) o IFLUW 
PUELETE GOL .3225+S0Q77CY1 2.1606 
C 4 ler 1 LAL 194 20/ els 22/ 234 24/4 PS 24s 
Kole} 7001800190 0200 02102365 e230 eee e250 e2hN) TFL IW 
*INSERT GUL .763 


NAMELIST /JITTER/S JITANGeJITOIS 
JITANG = THE ANGLE UF JITTER (IN MICRORADTANS) 


JITOIS 2 THE OISTANCE PRUPAGATED SUCH TraT THE JITTER 
S16mMa 1S JLITANG®UTTOIS®} .Eo& 


ANAN 
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: ae : a MT anes — ato snaitna : ; _— et 


c 
*OELE 
C oe 
Cc 

230 


231 
1336 


223 


235 
240 
250 


2en 


PINGSE 


@lse 


PINSE 


L746 
PUrELE 


TE GUL.327 


SHSSHSSSSHSSSSSHSSHSHSSHSSHSHSSSHSHHEHHHSHSHESHOSHSSSESHHHSOHHOHSHSHS 


[UTR 2 [yTQe] 
TE CANOT.LINIT) SO TO 231 

READ CINeUITTER?) 

ANC (SefUTRel) 2 JITANG®!].E-6 

ARC (GeTUTQIe2) = YITOIS 

SIGXY SABC (6 eTUTRel) PABC (Oe LUTR ee) 

WRITE (601836) ANC (GelLUTRel) eSIGAY 

FORMAT (45H 998@ SEAM JITTER MUDEL CALLED *¢*%e9TD OF VIAs 
XZ3HTION ANGLE (RANTANS) 2 eGleedeBXe%: STM DEVIATION (SIGAY) 2% 
XG12.5) 

MO 233 IJ21eNOA 

US(IJ) & CUR C2®I Jol) ee2 © CUR(2¢I 4) 202 

NX BR(2) © KCL) 

CALL JITRAGIDXeSIGXY) 

OV 235) [Ja] oNOd 

TJIU & 2eT) 

CuR(TuTvel) s5ORT(US(IU)) 

CUR(TUTY) = 0.0 

CONTINUE 

CONTINUE 


Cost Te 

fwuNAL -s 1 

G9 TA Gay 
ml GOL. 32 

FITVALENCE COSCE) oCEIL CLES o (CUMC LY eCuUCL)) 
NIMENSTON IS (146 3R4) »CUR (S27? Ha? 
aT Glels 

® eUSeCUr 

REAL JITTANGeJITUIS 
RT CYCLED.A? 

NO Lie TLeleNPTS 

eC{T) 2 XSAVE CTT? 
TF ‘WAL. 

SUBROUTINE QUAL (TRPHASE oS LSAVE IPL Fe TITLE oP e ANS eRe RFE e SI GANG) 


PINSERT QUAL.1] 


PULLE 
a3 
40 


late 


48 
79 


VATA PT736.1461992N6/ 

Th IAL. 1LO7 

LF USTuAnGe] F295 71070 945 
SIGAY2F OS [GaNnu 

whl TF (He1236) SIGANY 
FORMAT (45 800* HEAM JITTER MUVEL CALLED ®eeeeSTN DEVI Ae 
RKIGMTTUN ANGLE 2 ¢ G2U05 ) 
CALL JT THAG(NKSAVE 6 SIORY) 
tiMARBOS 

OU 64 Yale PTS 

Jl atlue) ONPTS 

YO AM [RB] tlBTS 

T72leut 

TF (Umar. GF .US(I2)9G0 TO Sa 
YMAKBUS (17) 

XPEAK ax (T) 

YPR OK ar (5) 

TAI Lsee!7 

CMT fel) 280K Ph sci Z)) 
CUR(TZ17) = 04n 

CUNT TNte 

MER BIIMER/S LON. 


INSERT MAIN, 300 


i 


SUBROUTINE GLTONG(OXeSTUAY?) 
TeIS SUMROUITINE MOOLTFIES Tee FaR Feu INTENSITY OLSTRISUTTUN 
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Fa 


ANANA 


2 Tatas a 65 
my ya 


MIOFL THE FRFECTS OF eat JETTEWS THE GT TTER TS ASSUMMEN TO 
GAUSSTAN, SINCE THE Rb TING INTENSTTY IS THE CONVOLUTION OF Tue 
UNO ESTURREY PUTENSTTY wil Tre GAUSS[ANs TmE QPeERATION IS PERFORME 
Gy THe FET ON FAaCH FUNCT [Ow ALUNEs® MULTIPLYING [HE KFSULTS. 

PEM watts Tee INVERSE FET. OVE eh/26/7H, 

LEVEL PeCiueCl 

COMPLE A Cuedl 

COMMON PMELT/SCHUL TSIM) 09133024) oK 6128) owl eNPTS ONPY CDRX DRY 

CIMMARS SCG/EL V7 993 

DEMENSTON Nw) 

NATTA BT 42,1415 937 

SPITS 2 NPTS @ NMTS 

Yea st 7a 


Malet) 


CI(M) = CMPLX (US(M) 9Ueu) 
10 PPw = PRWweus(M) 
NNN C1) sNeTS 
NND (2) sNPTS 
NAR = 2°®NPTSQ 
NPOZENPTS/2 
CALL FOURT (CleNARetNUe]? 
SIGEXP = 2,9(SIGKY ® PI) ee, 
SINEZ(X(INPTS) @X(1))/e0 © UAS ee 
OF S.5/5STDE 
AO 20 Jal enPfs 
YSWS(( Y=] ) eNF) e020 
TF (JeGT NPV?) YS 92 ( (UPNP TQ <1) MDF) 2e0 
IRS (Je) ONRPTS 
90 20 yetenPTS 
RSQ2( (Tal) eQF) aap 
{IF (T.GT.NPU2) xXSUS((L=NPTS<1) O0F) ooo 
Kalen 
ON CLURI SCL (K) OF AP (SI GEAR ® (ASQ 0 1SQ) 
CALL FHURT (CIleNARsNNUe ©)) 
Md 30 KK z) eNPTSU 
WS OKK) SCABS (CI (KK) ) /NETSu 
30 PRPWN = PPWNeUS (KK) 
PWRFAC = PPW/PPWN 
DU 40 MMe] NPTSY 
40 US(MM) = US (MM) @RWwREAC 
WRITE (e100) PWHEAC 
10G FORMAT (//5Xe8THE POWER HAS BEEN SCALEN BY A FACTOR UF @sGl2eSe 
X SIN SURROUTING UTTRKG.9//) 
ReETUAN 
ENO 


®UELETE MAIN. 60 


NAMELIST/QLOT/TITLE LOL T eOne ISAV es IPHASE pRAR ORE eSTGANG 


PNELETE MAIN,230 


210 CALL tataL (CTPHASE eISAVe TUL Toe TITLE oRER AS eNB RF eSTGANG) 


*INSERT “MAIN. 22 


NATA STGANG /0.0/ 


APPENDIX B 
SOQ USERS GUIDE UPDATES 
JUNE 1978 TO JANUARY 1979 


INTRODUCTION 


This appendix documents those changes made to the initial SOQ code between 
June 1978 and January 1979. The changes incorporated in the code are those that 
have become generally useful for the physical optics simulation problems which 
have been solved using the SOQ code. The users guide updates are also prepared 
to clarify and correct the initial description of the SOQ code as documented 
and delivered to AFWL on 1 March 1978, in the Preliminary SOQ Users Guide. This 
document supercedes previous written material on SOQ code documentation. The 


organization of the SOQ Users Guide Updates is as follows: 


Section BI New Subroutines 


2. Subroutine CPUTIM 
Subroutine LISTER 


us 


Section BIT Code Changes/Correction 


1. Subroutine ZERN 
BI. NEW SUBROUTINES | 
| 


1. SUBROUTINE ZERN 


Zernike polynomial terms give the SOQ code the ability to model mirrors 
with arbitrary surfaces. This subroutine also provides the determination of 


sensitivity of a given system to the level of these Zernike terms. 


a. Relevant formalism -- The Zernike Polynomials are an orthogonal 
set of polynomials used to describe phase front aberrations. The low order i 
terms of this set correspond to the low order Gauss-Seidel aberrations, such 


as piston, tilt, defocus, astigmatism, coma, and clover. A list of these 


polynomials, Z(k), is given in Table B-1. 


¥ 
¥ 


TABLE B1. ZERNIKE POLYNOMIALS 


1.0 3 | (4ar* - 3r*) sin2e Le 
Rcos@ R*cos4e 


Rsiné R*sinde 


2 


ar? - 1 (10R> - 12R° 


+ 3R) cosé 


R7cos20 (1OR” - 12R° + 3R) sine 


R?sin2e (SR° - 4R°) cos3é 


(3R> - 2R) cose (SR° - 4R°) sin3e 


(3R° - 2R) sine R°cos5@ 

R°cos36 R°sinse 

x 6 4 2 

R°sin3e 20r° - 3orn* + 12R° - 1 

on” - oR? + 1 7or® - 140r° + seor® - 210R* 
car* ~ 3r2) cos2e 252R29 - 630r® + seor® - 210R4 


+ 30R? - 1 


The phase applied is 


24 
Ag > 2nP 24 (R, 9) 
k=1 
= Ao(1,9) 
- = R<l 


° 


Rol (B1) 


If the Zernike radius Ry is specified to be zero, it is a flag to set the 


phase identically equal to zero. | 


b. Fortran formalism -- Subroutine ZERN is called by GDL with IFLOW 
= 24. Namelist ZERNS contains the Zernike radius Ry as well as the coeffi- 
cients of the Zernike polynomials to be applied P(I) [I = (1,24). 
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ANNI >, aptnene cba pase SPF ; mee oa » iat 


Due to excessive use of the FRINGE program, one can also input fringe 
coefficients (PFRNG(1)), corresponding to the 24 Zernike polynomials to be 
applied. The PFRNG coefficients are converted to P coefficients in subroutine 
GDL. 


NAMELIST /ZERNS/ RO, P, PFRNG 


Argument List RO, P 
Commons /MELT/ 
Externals None 


IDENT ZRNIKE computer printouts follow. 


Le iT Peay 


OTe st she En 
toe 
wehbe Te VETTE mets 
an WAS ene ee ore ee) 
mere Té sltte~ eer 
fat z Ys Awest vote The Aa Few bX ES Poe iTS OF WAVES 
hime a THR TT Te a oe 
PENSrOT ITE Th wean 
LMT AL Bap rreeb geek Vt bere A oer Tem td IISA RIA IN 
OESSE QT GETTER pu 


TOPOS SEHHESSEHES HEH EEHEEFESH SHES ESSE HEH OH HOSE EEHEHHELESEH ESOS OHSS 


c AMOL y Cb hae 
Cem HE HROH OHHH ESE HEEEHE HHH SHES HERES SEHHES HHH HSH ESHEHESHEEELE SEH EHOO EY 
Aor [shes Bz T/tw- & | 


th Cowl et ff) wy F> soe 
Fue { eK 2 ,F Ny 
: roy Pad Tepes 

Pwr PLT) Bory 
WY Jud [te dn 

Dow MED ast T} 2 oy 
WE ANY Ce 7h IG?) 
ty 2a Taledias 

Prd FR CPr mts rl Pre Peded Fal veh ae ly, 
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SUBROUTINE CPUTIM 
Subroutine CPUTIM has been activated for the CDC computer to print out 
the amount of CPU seconds used by the kinetics package, which is driven by 
Subroutine REGAIN. On the Cyber 176 a system routine 
A = Second (B) 


returns the CP time since start of job, in seconds, to both A and B. 
FORTRAN: 
Argument List: 


IT = 100* time since start of program 
Commons None 


Externals None 


IDENT CPUTIM computer printout follows. 
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3. SUBROUTINE LISTER 


~ 


Subroutine LISTER was activated so that the output of the resonator. 
design program RESDES or an arbitrary file may be read internally and 
reprinted in the output of the SOQ code. LISTER reads an 80~column file, 
designated as Tape K, and reproduces it in the SOQ-designated system output 
file with pagination defined the same as on Tape K. 

FORTRAN: 
LISTER is called anytime IRSDS, is nonzero in namelist START. 
Argument List: 
K (= IRSDS from START) 
= tape number of the file to be replicated 


Commons: None 


Externals: None 


IDENT LISTER computer printout follows. 
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15 RPE«wINf K 
ar [ TE (0 aS) 

VW Fw AAT CTE LD 
HE TH Iiiy 
emf) 

“ar 
OCISERT MSTNA LSS 

TF CT8S0S.NEeU) CALL LISTFERCIRSIS) 
IF (awl el Fete) AEWINM 50 
T2SnS5 = 19 


BII. CODE CHANGES/CORRECTIONS 


The code modifications and corrections included in the code are 
described below by their update file name. The reason for the change, the 
structure, and the listing are included below: 


1.  *ID SOQMAP 


This update provides a cross-reference map to the SO0Q79128 code. The 
first section lists each routine in the order of appearance in the SOQ code 
with its commons and externals. Also given is a list of all routines that 
call it. The second section lists every common block in the SOQ code with 
the subroutines possessing that common block. 
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ra *ID ABCMAP 
Current allocations of the ABC (I, J, K) array are presented here for 
ease in future updating. 
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3. *ID PLTFIX 


Ident PLTFIX modifies the printer-plotting package in the SOQ code. 
This new plot package: 


a. Prints DCALC, IMAX, DCALC FLUX along with the location of the 
center of the beam (DRX, DRY) and the bottom of every iso-inten- 


sity plot 
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Prints a blank for every value of intensity less than 0.01 
*UMAX (UMAX is maximum intensity) and puts a border around 
the outside in column 1 of NPTS and row 1 to NPY 


Allows for selective plotting, based on the new namelist 
parameter KPLOT in namelists PLOT and QLOT. 
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4. *ID ADDPRNT 


$ This section of updates was included to add information on intermediate 
a printout to CAVITY, STEP, GDL, and TILT: 


a. CAVITY - The incoming and outgoing total flux at 


each gain/phase section is now printed. 


b. STEP - At the beginning of STEP, current values for 
DRX, DRY, RAPTR, NREG, and WNOW are printed, and the 
incoming flux calculated. At the end of STEP, mod- 
ified values of DRX, DRY, NREG, and WNOW are printed 
along with the percent flux lost during the propa- 
gation step. This last parameter (percent flux lost) 
indicates how much of the beam has been propagated 
out of the calculation mesh and, therefore, lost by 
windowing in S-space and K-space (Fourier Transform 
Space). 

c. GDL - At the end of any IFLOW the code now prints 
out total DCALC FLUX, DCALC, and the location and 
magnitude of IMAX. 


d. TILT - Subroutine TILT now prints out the mirror 
radius of curvature necessary to remove the beam 
radius of curvature found by TILT. 
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Rw shies ® (X12) 0K (E)) SOPPFLOATINGTS/NPYI SINUS 
LF (CUE Ee Le weit Fie) POWS = PO R/WNOW SO? 
@T ISAT STEM. 2% 
Or 2 NRT S@ RY 
i) 40) Ta) oN 
tf = 2] 
oop POeS = Pied © CIO(T Te), ee? © Cur ( TI) *#2 
BPNwA = PNWAXLK (2) OK 1d Eee POELOATONMTOSNPYISINONN 
LF ONME eh ive Lee NMF SP oe 7) PNWA SF POWA/SWHIWee? 
DELLE Ss (BO WASPDyA)ZONWreI AN, 
LE CORT I Ee De Meet IWE GE oD) ANI eo TPRNT otk oth) WRITE (Ae 3000) 
XK AXOF AL ce IVVRE AL MRE Te WN IWe DELP 
PINN FURAATI SGX eo ME RTT ING SURAMITINE STEY. CURRENT PROPAGATION®, 
4 & COBWwaMETFERS: 


X JH 9 OCW z®eslrebe 

: XR /HX,eONY Beal eeee 
g SRL, ONBRE see THe 
C/E, owhOW BP oelrieewe 


Be S4X,OREQCENT FLICK LOST s®%ebl7e4) 
OELETE STrP.1L91 
Te FURMAT(/aX 0} 9H STHE MM TNTENSITY 26Giec eS) 
@INGFQT STEP,2A5S5 
wNOW 2 1.0 
LE CTORNT NE 91) SRITE (He WINN) NERF AL oe YREAL eNHF Ge wNOW NEL 
OTHGEWT STEO AKG 
wnQOe = 1,0 


GO 


ONELETE GL eHOT eG IL AGA 
uU4ax 2 U,0 
KMAK @ A(]) 
. YMAX = X()) 
DQ 7R Yelenry 
JL 3 (Ue) ) ONRTS 
OO 7R Psy esye#>rs 
1Z = Tel 
AYINT & CICT) OCONUEICUET 
IF CUMAX.GT.LAYINT) 69 TO 7A 
MAX 2 XYVINT 
AMAR 2 X(T) 
YAAX @ KY) 
74 PPw =z PRWeXYIN’ 
LE INWE GH eEQ SL ON RES Fie) UMAXBUIMAA/WNO WOO? 
UIMAAK @& 1iMACs10049, 
= WANMAX B SORT (AMALPODOYM AK OO?) 
tan PELE FE GIN. .AOAR 
TF (MSTEP NF ol) Sat TE (079) PRWK eNCALUR oe UMAXK o XMAX & YMAL eRADMAX 
®NELETF ol 870 , 
MUX BweGleeasArne L2HICALC SoFHe2/Hae lLAn[ Max SoGlreselUXe 
M2NHLOCBTED AL (CAeY) & Celts LeoGdeltee's] Powell) s 
XY SGAPKeGHAT RAN[US eO127.9) 
MYELFTE Gil 47} 
IF CMSTER CEN. 1) WRITE (Ae TTY} HO MK eNCALCPetIMAXH eo XMAX eo YMAX eRalimMar 
*HFLETE GOL A723 
MUN Bets] 2 eG/Ane LAHNCALC SeoFRA/HA eo LOHIMAX BeGleewel My 
MPNMLOCETEN AT Cee) & CelsLPeselteelslcoSelmde 
X SGAX GMAT FaANTUSeGLAes) 
PNP LETE GOL eS1AeGin o51 4 
{QN4au = 1 
60 TO 999 


TIL? 
PNFLFTE CYCLE 9.293 
TauRPl = 2, @ranciiw 
LEC TPS .Gr 7) aT TE (4667) RAC IR e TWORLC 
OOFLETR CYCLE9.295 
KLAR eo SAMMHAGE FOONT CHIRVATUIME 2 RANICIIN BeGlLPA eee te CM/ 

JY%Re® NOTE © Tes CumvaTiiRe CAN AF REMOVED wIfH A MITRRORe, 
PVG oF NOT WANG welt eee IIMS ZeMmaANCUR AS DEFINES) GROVE © / 
10x 6@ SEGATI VE VANCUMW TS A CONVERGING PHASE FRONT del CHOes 
VPOne® CAN 6 WEMyEI VET A CONVEX (NEGATIVE RANG) “[YPRUR 66s 


an Ma eA 


x 41) 


Ca mE ae oe 


5. *ID SCLPWK 


Ident SCLPWR modifies the IFLOW = 12 section of GDL to allow for 
scaling of the beam to a specific power TRANS. 


[FDRNT SCL PW 


OT PUT SCL Pak 
HM, 
PLUISFOT FOL atm 

MAG @ 1, 
PELE TE GOL eG WE 249] 
© eee TETAS .LE elo) THE FLFELA TS SCALE! BY SaRTCPRANS) 22 aG 
Cf Pee TE CTOANS GT] 60) THE FIELD TS SCALED TO THE PuwER "TRANS 
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| eee 


I] PUL = OIC Che MET e a) 
AMals 3 BC (A etLT ea) 
TRANS = BOLD 
FP CTOARS oth etet) GO Tu ays 
Wipwa Sf UA0 
UO FBSA L7s,eude 
We PNEw Ss RPNE Ww © CTL) OCOMIGICUCTZ)) 
PYFW FS PNE eh OCR CPP MK (CL) } OO?) O(NPTSSNMY) /S1000, 
PF UNWE (yo isg)] eV e IWEG ENG 2) PE WSONE W/WNUWOe? 
GAO LF CTREPS GLE Le} ONEw = AMAGOOD 
STEANSSSGRT (C2 IL I sOME Ww) 
weT Te (60357) BYC(CLoeMEL ToS) oAnC(AeMEL Soy) 
ST YSERT GN 2a94 
TF (TranS.G°.1) [4n4ta] 


6. *ID TBLUM 


Two errors in subroutine TBLUM are corrected by this ident. The 


following listing is self-explanatory. 


[Meer Tatas 


Orel Tet iim 
Taya 
SUELETE Pal yuma 
TF (MALO eGFete) ARITE (4 e595) AXTAL 
SORE TE TAL VM .anr 
CVT 2 (9A SAAS HAL FNS OMT YS (werd OCPOl) 1) He (1 6/50) 


7. *ID REMSPH 


Ident REMSPH allows the removal of defocus and/or tilt using a call 
to subroutine QUAL, and to continue with this optimized beam. This optimized 
field can be plotted and written to a local file specified by IWRITE. 


IWRITE. FT.0 sets IW = IWRITE 
IWRITE. LT.O sets in IW = -IWRITE and returns to SOQ 
immediately. 


If desired, the non-optimized field can be read in using ISAV = 1 in 
namelist QLOT. 


PF at May spry 


PTR YT WY SP 
SAL 
*#TSeef tac. 
(FCT aE Eee) HD TO Mic 
wide S FRI G@ypy 
To s feamt fe 
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TE CT wwRITE LTO) in = ml we tte 
we LTR CLad) (CrCl ape dase e Nie) ek ewX slr eN[ Te Saver 
GEwTNi Lo 
wi [TE (Reh d) T te TPHASE 

SY FOR FATS eke PCI) AAS REC, AWITTEN ON UNTT Mel 30% FROM QUAL Se 
X @ a] Tr [Pease sHeelPss) 

TE (CT aRITELGE 9) 8) To ad 

IPyuT = | 

TE (Ke ET GT et) wRITE (Ho WN) NETETE 

TE (KMPLOT.GT.2) CALL FOLOT(KPLOT) 

LE OTSAVE Rei) WR ONDE 7S (COC ND e EXS be) OX gINX CUNY 


AF alet) 7 
Wwe Trine | 
AO CONT Tvue H 
“aa [8 ig 
@QEL FTE MAING CP? “4 
200 KeOLOTSL OAD £ 


HEADS eQE OT) 


*ID CHGNPT 


Ident CHGNPT increases the flexibility of IFLOW = 6 in two ways: 


a. Reoverlap the beam, letting the code find the 
original DCALC by setting DIBEAM = 0 


Change the number of points in the beam, by interpolation, 
by specifying NEWNPT and NEWNPY. On a subsequent call to 
START, set NNPTS equal to the value of NEWNPT or NPTS will 
be reset to the previous value of NNPTS. 


a oe ee OTT Sa 


OTORET CHONRT 
AOL 
OUELLETTE (ote MY 
150 [Ler AastweA o 1 
TF Cew Teo Tl) tod TO Peg 
bm RVC I Neh Ct PV) 
ANCL Lah Aen) @ atyede? 
159 Wheel) 2 APC Loe TwEAeh) 
*T"SEWT CYCLES A 
NPTS 2 NF anor 
NPY 2 PbhwARY 
t NOM 8s MPV OeRe TS 
PPuSFOl CYCLE I.4 
KN ih TP ah eg [AL WSN 
ME IPT Ora SE dS YY AWE THE NE STRED NUMBER OF BOLYTS 
SA WE A Tae ELELO, 
fHfOyyw 2 Peto Cre Tee Cut, OF FAULTS TO LMIn = le 
T1 FS Use TO RENORMALIZE THE gAWwe Re SONWATUR FEENRACK FIELNA. 
UNSFOT Gi 2A? 
WENT 2 on 
Nery 2 fy 
IMIWGM = } 


“aie siren eaceahipetaet eas ons area be waelmstie dy 


i in 2 el see oh aps Sept 


eeovany 


aged 


PTNSFRT Gol ema) 
aC lmebel) s&s Vey 
AGC (Pe lol) 3 Ae ahwy 
ANC ladeleo}) 2 Tepes 
TF OVD PRAM FI.) MICK LeLel) SF CXCAT ORC TD) MEL OAT INPIS) 
Mie EE TE aN gar? 
Newel s AUC (Hele tlo.t) 
New ley =F AeCiPelelleoedil 
[412s 2 AAC (Asletioga) 
TF (NG ate TF OGD) ME Whe TSNRPTS 
TP ONE AMP YF) ME NPY SNRY 
ROFL Ss TM SsNE at Tad, 
POELE TE wt ene a 
a AA [ie Ae sh (NPT 
@NELE TRE GIL Ane 
NO AL MYR Leute sary 
SLE TE SOL HRA 1 | 
UM 4 MKa] ob wT 
PORLE TE AMVAAST? i 
NE width? oS oh yitP here Wd py 
AMNOWSQ S 1, 
LP CME ek dele Me elIWE Ie F 27) WMI ISQSwhUwee? 
PwA Sf ROWA/S aN WS i 
Yr Ss ywhit J yh! pwS ) | 
WY G24 Pes] ei aynu | 
{ 


9.  *ID MISCFX 
This ident corrects minor errors and adds two parameters to namelist 
START. 


. 


a. *D GDL 384,385 This change in format compacts 
this part of the printout to 80 columns for 4 5 
or fewer struts. 


b. *D GDL 884, 885 This change removes the S from 
column 1 so that output can be put on microfiche. 
It also corrects an error in the BARE updates so 
that the CU field is read in at the end of a con- 
verged iteration. 


c.  *D CORR 1.23, 24 This change removes $ from column 
1 in the output. 


d. *D JITTER .83, 86 This change updates the indices 
in the ABC array which were defined originally in 
reverse order. 


e. *D BARE .11 This change corrects the size of the 
loop from MUT to NOB. 


*D Cycle 9. 119, 120 Previously for IPS = 2, 
the iteration counter KOUNT was not updated. 


*I Cycle 9.99 Focal = 1.E50 defaults the radius of 
curvature of the beam to "infinity." 


*f STEP .40 This change activates the IIPS # 0 
option in STEP. Setting DELZ = 0 allows the 


removal of tilt and/or calculated sphere without 


propagating the field. 


*D GDL .827 This statement was redundant. 


*D BARE .86 The parameter RGAIN allows 
the option of not calling REGAIN at the end 


of an iteration. 


The parameter IFLGAP is included so that aperture loads are printed for all 


apertures in the optical train. 


the at ESCO 


*TE ER af “TSCA 


THIS PV Cm FK OT AP ee FOR G ALA ANNs MAR AMETERS Tr) 
vaVEL IST STAWT, 


Gf, 
ST LE TE lo Wwe sg SS 
TNO, IF SPR ITS ee LreF tal At KV CRANE Sets lie Selden] otal AXe 
PV Ag TA ARTE ee 1 oe he PHTHE TASS eAULN SY) 
POPE ETE GOL ieee ZO GURY 
FAD EN AATIs/ Be LEG 44) 77 33N ITEPATION [5S CONVEHGED « 
t mH NMP frivelLaelor [TRAATIANS FAVE e tp LA l4aey ss) 
ARAN CIP COUHCT ZA) oe PZB oititPR) eX el ArersiWY 
WRaTNd % 


LISTARO 
HELE TE C RR] AIC WH] 07% 
SEC eAHCAw eT ee LOC TD eo bir CPAP) oe POC bers) oe LUC LH) oe LUC LMS) « 
I SUV Vem) OLE LZ) ob tee S72 CIM ee OK OL LNML 2A Ja HG) 5K 


ar 


PORPLE TR TT Tew em te JP PPR eae 
ARC CL eLUTo eat = UP TANGFE] Fen 
AAC OAc] ITH es) = YTTITS 
PAL STety SAMO OLe LUT eA) MARE (Ae fl IT em) 
wep hr (mel tay ANCL ae LITK em) ea LGry 


Cavity 
ORL RTE WAawF .)] ) 
139 wATTe (7) 


see qi rT AE 


COHOT2Z) aL 7st adh Or} 


rict 
OUFLETE CYCLE MeL Tee CYCLE 4.120 
eS <a st eto] 
TR CPOES Ewe) HT) Se 
PL USkaT CrCl ba ee 


CCA = Pek Ss. 
STFS 
PP uSFOl STEP, ou : 
we ANClin S €2 10 Aw ae 
CF CDEL Zeb tee) RET Sey ‘ 
SO 
PPELETE BMT 
POELETE RAWF RA 
| IF O98GR TINY CAaLt SEGATSN(SCT ITED) 
“At 
. OIL VSREST malwel 4 
_ LOGICAL #GATN 
= PINSFRT AAR OY 
| c T2S98 [5 Tae TARE “HidReR OF TRE AN@CULUMN FILE TO BE COPTED Ta 
) is ceiTw IP Ay LISTE, fF [eSnSsue0e LISTER [5 “OT CALLED, 
c RGOAT™ = oh alLSFe TURNS OFF THE CALL TO wEGaAIN IN TFLOWws7, 
PLUSFYT 4A Ne} Se 
Wisagth s eTwur, 
TRS05 = 4 
IFLGAR = 4 
SINSR OT ANINS EAL 
TFL sar = 3 
SEL ETE adwh 5 
Ro LAAWE PLOTS eRGAINe TRSNS 
OTUSEUT MATNS? 
CMAN SSVTY 4S RGAT Me TEL GAP 
AnL ‘4 
@T sSFoOr sore. l7 , 


COMMON SGVTYMS AGAIN TFLGOP 
PELE TE GOL. KH} 
TF CT CuTL cE de Le ANN STEEL GAM CFI.) GO TO 96H 


10. *ID FXQUAL 


The quality program has been updated to include more options and more 
printouts. See also *ID PROP and *ID RMVSPH for other additions to QUAL 


a. IPRNT. This parameter was added to suppress the 
additional STEP output (from *ID ADDPRNT) when STEP 
! is called from subroutine QUAL. It was also added 
to namelist PROPGT for the same purpose. 


b. The output of the focal plane search was modified to 
: print out more information. 


Ci Additions to QLOT 


RBB (New meaning) 
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IWRITE 


PROP 
IRYFF 


KPLOT 


I TABLE 


ICTRD 


(1) 


(see *ID RMVSPHP) 


(see *ID PROP) 


RBB: 


If RBB is input as other than one, 
QUAL will find the quality infor- 
mation for RBB (RA/D) 


IRYTFF.GT.O writes far field 
to unit IRYTFF 


KPLOT.GT.O plots the far field by 
calling IPLOT (KPLOT) 


ITABLE = 0 finds quality table and 
plots information 
= 1 does not do the above 


ICTRD is used for ITABLE = 0, 

= 0 chooses the optimal focal length 
based on the highest 1.0 RA/D quality 
about IMAX, then constructs the quality 
table based on the better of the two 
beam qualities at that focal length 
(Default and same as previous). 

= 1 calculates quality table about cen- 
troid for optimum. 


= 2 calculates quality table about IMAX 
for optimum. 


= 3 finds the optimum value about either 
centroid or IMAX chosen for the highest 
1.0RA/D quality. 


poet F au AL 


@TOENT FROQUAL 
"uae 
@TuSFQT JTTTew. te 
C ete MRE TS Tee FAW FEED CRNUTERLINE INTENSITY DUE TO A PLANE wave 


c APR ATUME TO A TAME TER em wht a CONVERGING LENS UF FOCAL LENGTH 
Cc F ABULTED af Tee WE ae FIFLA, Toe PT al, PowrR IN THE APEQTUREN 
cr PLANE wAVFE ITS Tee Sa4F AS THAT OF THe CURRENS CU FIELD. 


PENSFOT CYCLE Y.17 
COMMON SSTHOWI SZ TRANT 


BML 
PINSEWT GOLelee 
x 9 | PRNT 

@INSFer GOL,138 
c 
¢ [2eNT TS A FLAG FOW PRINE TYG NREG AND wNOW FROM STFP 
c 2 0 AONTT OBRINT | 
Cc 2 1 PRINT (OFF AULT) 
@TSEGT wt esl ! 

[PWNT 3 } | 


PINSFAT GOL .S9N 
BBC (He ISTE Ss 3) 3 TRRNT 
@INSFRT GOL 2593 


{ 
LIPRWE 3 ARC(Ae1STER63) © 2901 { 
PINSEST (bela 
COMMON /STRPQwI/s TPORT 
eek 
®TISFRT JITTERY. LAA H 
TE CKPLOT SOT. 7) WRITE (6021009) FFTIFL 
WAN FARMAT (LHL een Aass) : 
TP CKOLUT NE.) CALL [PLOT (XPLOT) 
[PeNT = ] 
IF CTRYTEE Nb ot) weITECIRYTFF) (CUCTK) oe IX2L eNOR) OX eORXODRYONIT e SAVE 
LF CTRYTFF ONE.) APYINU [RYTFF 
TE CTHWYTFE NE.) WRITE (AAD) [RYTFF 
MAN FAWMAT (LOX SPR AR FIELD AAS SEEN WRITTEN TO UNTT@e] a) 
STFR 
STNSERT STEMS 
COMMON /STH wi 7 TRPRNT 
Cavity 
eT vSeet CaviTy.9 
COMMON sSTRPUuwl/ TRPYNT 
®RNSERT CAVITY.102 
[Pwiy 2 } 
Qa 
PINSEMT JITTER LIIG® 
HATA FRTITL 41 4%64 eGHE AM e4HF TEL conf PLesnOTS «2%6H 4 
Data NETITL /14%@4 eGR TT GHAl ZF eahi) FL eGHFLI) 0 P Pam 7 
ATA SAVE 4199067 
[PRNT 2 0 
wAl = 1." 
wis = 5,7 
@PINSFOT Wwthdle leo 
C @ee [SAVE = 9: PENT) [ta FAM FLELO FROM tie [T a, 
@LUSFOT wields 
fc eee [Save 2 1 ot SAVE NF aw FIELD ON UNIT 7. 


@TNSEIT QUAL 2] 

C #ee [SAVE ee] 3s WEAN EAR FTELO FROM UNIT Ye 

@PUSERT CYCLES L2n 

C oe wRITE CH FTEL'Y wETr LENS APPLIEI (FOC AL LENGTH F) TO UINTT be 
OnFLETE WAL ell? 


CALL PO aww (NETS cA eR elTe XPE AR 9 YPF AK CRE] p PHM) 
TE CISTERP eh eh) CALL POWdLW INP TS oD X ek eUS ce XPEAK » YRE ARK RNS ePRRS) 
[TF CISTE RP LFO RS CANT CRANE OL eo?) 
KR CALL PUwWANe (PTS eX eA CUS ce APE AK 0 YPF AR RH ePRRRA) 
PNELFTE QUAL LAL eQiAL.1A/ 
WRITE (Hel 37) Ot] Pa eXCINT eC YCINT ce RH] oe PRK clLIMXK gp XPFAK eC YPF AK 
x PWSAVK DR eS TRF HL 
TF (Re .NE ele) 
KR AW TE (el 37) BH ePRAMHM XC INT Ae YCINT eH oe PRRRGK oe UMAK oO XPE AK Ce YME AK» 
x VANSAVEK NBeSTRERL 
19 FORMAT OS/156 SCALC FLUKA IN eF5 e220 fH wL/ISeGlAede® AROUT CENTRO, 
KLAR eLUHCONENTYATE Se AGL 2ee/1 SM DCALC FLUX IN F502 06H RL/NZeG1LAede 
XLOH AAOUT TMOd DF ed 2emel 24 CNORNINATES ce 2G612.4/13H TOTAL DCALC « 
XSAPLUX BORLA cheb ke 72h WEFERENCE NTAMF TERS eF GH e20/ 
X }94 STHR tL TMTENSETY 20011 04/) 
ANS (WT 2 PRHSK/PwSAVK@1LO0. 
BOSCNT 2 PQASK s/PWSAVA e109, 
WRITE (06010) RQOSINT eHUSCNT 
SOLO FOMMAT(SLOXSANOTES CENTINEIN an [Max COORDINATES arf ING, 
© @ CENTIMETEHSe@ss/©@ Q BADUT [MAX FOR SRL/UBOeGl2e40 
2 10K BY ABMIT TRE CENTROTH) FORM SH /N290612.4/) 
PNELETE tual el12 
CALL POwwOW (HOTS MK eX eS eKCING oYCINT eRH] PH) 
ZLNSQ = ZENe7tYD 
PINSENT otalLe lo 
IF (TSTEP.NF mH) SH) THY ANNA 
C eee FIND Bier Ty aR). 
CALL Pl awOW (PTS IX 0K othSe MPEAK Ce YPF AK eRAS ePRRS) 
CALL PURWID TINPTS oe Xe keUSe MCIN oO YCI NT WAS se PRA) 
PHHS = BPLNGE/L1NS) 
PORSR = PRIS/1INAN. 
PHES 3s PWHSOe7L DS) 
PHRSK = PARSsSLYNN, 
TE (Re OeD 6) GO TO Anu 
FIND PuWwFR Iw Q8eRL ZI) 
CALL PUWWOW (NETS eJX ok otiSe XRF AK oe YPE AK eRe PRANK) 
CALL POWWI AI (NPTS So X eS Oo KCINI eo YCINT eB ePRHWA) 
Play S OVQQNeZ1 NS) 
PORQHK 2 S9ORHASLOANO, 
PQAHYR 2 PUNWOAP71 95.) 
POHPHK 2 PONVA/sINOAN, 
C eee WETHION TO CENTIMETERS FOR OUTHUT, 
2900 KRCINT = ACTNTOeZLA 
YCInT 2 YCINT@eZLED 
APEAR = XPFAR®ZLN 
YPE AR 2 YPFAK®7LN 
PHIISTEP) 2 PAK 
HQT 2 PRAJPYSAVE*LNU. 
HOQOCNT & PAR/VNSAVE LN. 
LF (UROP NE 69.9) GH TO 340 


POLE TE CYCIL& 4 eGae CYCLE 7040 
TE CISTRR EDL) we TR (Ae 59] U0) 

Sap FP ywAar (S73 466 30M FLuxtet) IN PORL si ABOUT of 
RC LA TRIAL FUCA oe AU eSHPOTAL enke dH TMAX cENT 
X AX eae ShwF ty, af 
¥ {ee LEASTHS etc ese FL ikem Ke TOM (APF AK oe YPF AK) (XCTINTeYCINT) 
NM Ae VINE ST T Ys 1 Re IG ( ]4e)) : 
WRITE (05970) [STFU SE ePWSAVK oe PAK CANINT ePHK OHOCNT ep STREHL « 
« APE AK SYRE AK eXCINT A YCINI 

Sap POMMALT (SAW FellelaweGlAcGeAnehTeAelaal Hy eFTecelHloFanlel4) oAke 

CEP tlePae] eA) eth eF2,Gs JIKe lH Cob Me Jet Heer Hoteles 
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Y LKeL CePA, de ltta oF Metelh)) 
orFLETFE CYCLE9.17 
Deb HM (4) oh (09) PRIA) OXSAVE (127K) eFETITL (20) eNFTITL (20) eSave (10) 
PNELETR CYCLE DAT eCYCLE 2045 
PORTH = @190, 
C @ee FYI LUCATIOG OF MART MEM QUALITY AKROUT [MAX. 
NA V79 Lele 
TE (PCL) ott ePrOr) 39 FO 370 
POPPF s PIT) 
sv = 1] 
Wa CONTINGE 
C tee FIND LOCATION OF 4a4TMUM QUAL TTY ARDUT TRE CENTROTDe 
Mm) 375 [ales 
TF (PRCT) LE eRe) 6) TO 37S 
POPTH = PRT) 
{SvH = | 
275 CONT TNue 
Cente NETFRMENG FICAL LENGTH FOW OPTIMAL CALCULATION 
IF CI CTH ch eM eWe (CTR CE?) TOPTsISv 
TE CICTRE SFO el) TAMTSISv4y 
TF CUCT RI eRe 4d 30) TO JRO 
TORT = TSv 
TF (ROOTH GT PORT) (9PT = [SVb 
IAN FORTS ham (CT OPT) 
PNELFTE witalel it 
LECT CTR cE Met oe ( LICTIN cE Ve Vee ICT HH cE 0G) AND ePRALGT.PRA) ) 
x 50 TH 49 
ONFLETE tal e136 
5a TF (TTantL& Fuel) 30 TO 345 
WRITE (Heo) ACT NWT eo YCINT oF 
SS FORAAT (S20 0@ THE FOLLOWING QUALITY ToRALe [5S FOUND AROUTS, 
Ce PYIMITMATES (HeGLAebvlteeGleese%) FON F Sots 2.4) 
XCINT = XCINTZALEYO 
YCIN?T 2 YCINT/ZL9 
CALL PLTOT INP TS cf K eX UMA oe Se eltISe [PL To 
SOELETE STTTEN.I1M} 
Suazo rime Qral CLPHASE e [SAVE e LPL Te TETL E eRH eC ANS oDReRF eo SI GANG se PROP 
x <PLOTe I wHITEsLTAnL& se ICTR IRYTFFeNIT) 


MALS 


PNELETE JITTER? 
MAMELIST/QLOT/TI TLE s [QL T MA e LSAV es LPMOSE oKAR ORF OST GANG PROP > 
£ MPLOTe [wRITE eT TAYLE ICTR THYTEF 

@INGPYT MATN.HO 


ec PRIIO 0. CPRRFORMS FOCAL LENGTH OWTIMEZAT LOM 

c PROP LET .A. CALCINATRES Qual I Fy FOR TRE NOMINAL FOCAL LENGTH ONLY, 
r PROP LL Tet. CALCULATES QUALITY FIR TRE CHOSEN FUCAL LENGTH 

c (Fog ©PQUP) ONLY. 

c Two l TR etl oO SETS [4 2 TwRITE 

c TwRITFE LT. SETS 4a aelawiTF ANt) WETURNS TO SOO IMMEDIATELY. 

c TRYTEF 66T.0 wef TPS THE Fay FIFEL') fT UNIT [RYTEF 

c RPLOT eG Tea wMEOTS Tat baw FIELY AY CuoLLING [PLOT (KPLOT) 

Cc TPFaaee 2 8 Frys WALT TY TARLE af PLOTS TNFORMAT TON 

C 2 1 DOFS NOT DO TRE aROve 

C Icihany 1S JSF) FIR LTARLE = Oe 
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11. *ID FILAPR 


Ident FILAPR increases the generality of the aperture routine APRTR by 
adding filleted apertures. 


The outer fillet works by putting a circle in each 


of the four corners of radius Re (input as RFOUT 
or RFMOUT). 


The lightly shaded region is removed by a regular 
rectangular aperture, while the heavily-shaded por- 


tion indicates the region removed by the fillet. 


A central obscuration can be applied in a similar 


fashion. The result is a rectangular aperture with 


rounded corners: 


The input was a list of names for RFIN for 
namelist APTUR and RFMIN for namelist MIROR. 
P Subroutine APRTR has also been modified such that it now prints 


out maximum intensity on both the central obscuration, as well as on 


outer aperture. 


The bare resonator normalization aperture has been generalized to 


include the fillet as well as being any particular mirror number IM. 
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12. *ID NUDISKT 


Ident NUDISKT modifies the two I/O IFLOWS in GDL, IFLOW = 10 and 
IFLOW = 16. 


a. IFLOW 2 10 Two new options have been added to this 
IFLOW. Multiple fields can now be written to the 
same file by not rewinding the file between writes 
(RWIND = .F.). A file can also be written that can 
read at a terminal (READS = .T.). For this can the 
file is written in the following order: 


TITLE, NPTS, NPY 
(X, [x], I = 1, NPTS) 
DO 141, J 2 1, NPTS 


141 WRITE (IWRITE) ccuft + (J-1) vers], ta 
1, NPTS) 


Symmetric fields are unfolded before being written 


to tape for READS = .T. 


b. IFLOW * 16 This IFLOW has been updated so that 
formatted data can be read in as well as written 
out. The format has been modified to include more 


digits. 
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13. *ID CY4KIN 


The capabilities of the numerical SOQ kinetics package have been 
expanded to include oxygen, hydrogen, and R-branch transitions (9.4 uy oand). 


a. Oxygen has been upgraded from a structureless mole- 
cule to one that has structure. Therefore, there 
are now kinetics rate equations for the interaction 


of oxygen with the rest of the molecules from the 


b. Hydrogen has been included as a structureless 


{ 
| combustion process. 
} collision partner. 


| 


Cc. Previously the code has used P-Branch transitions 
(10.4 uw). Using GFACT less than 1 now activates the 
9.4 u R-Branch transition. 


In addition to the above major changes, two small additions have been made: 


(1) Input the Gladstone-Dale constant GDC in name- 


list CAVTY2. 
‘4 (2) Account for the gain length by the factor ZFACT, 
i also in CAVTY2. 
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14, *MIRFIX 


The MIRROR subroutine has been modified to calculate the effect of 


power-induced surface curvature when mirror reflectivities other than 
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the design value are used. The parameter, 5, is modified to change the 
center to edge distortion as a function of the mirror reflectivity. The 


parameter, RFLFAC, is used to scale 6 as input through the relation: 


ersten) (s) 


R = Mirror reflectivity 
P = Incident energy 
d = Design value 


Further, the MIRROR routine has been updated to include the calculation of 
its own value of mirror flux-induced distortion factor when mirrors are 
encountered off axis as noted by PHIAST # 0. This update has not been 
activated, since it would mean input file changes for all users. It is 
included in the code and will be activated by each user, when so desired. 
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15. *ID PROP 


The SOQ code calculation of far field performance is based on the 
analytical equivalence between the Fraunhofer pattern and the propaga- 


tion of a distribution with field curvature, f, a distance Z = f, using the 

Raleigh-Sommerfield formulation of the diffraction integral in the Fresnel 

degree of approximation. The SQQ far field calculation propagates the wave 

distribution, CU, a distance f, determined in a manner which preserves the 

correspondence between near field and far field coordinates, while accurately 
e resolving the energy spectrum in far field coordinates. 


In certain cases, however, it has become necessary to propagate the 
distribution CU to an arbitrary focal plane Z, using the SOQ calculational 
procedure, in order to obtain the effects of beam jitter at a fixed distance 
Z and to obtain the far field information scaled to same focal length. Since 
far field calculations are based on the use of "vacuum" propagation, the far 
field at any plane Z" is simply the scaled distribution at any other plane 
Z'. This can be shown by comparing the far field distributions in terms of 

the Fresnel integrals at two arbitrary focal planes Z' and 2'', where a field 
; curvature of f' = Z' and f" = 2" has been applied to obtain the distribution. ' 
| Comparison of these two distributions for the same transmitting aperture size 

leads to the following scaling. ; 


2 z 
~ik (f - 2) ik +7 (1 - 3) 3 
CUpw = CUzr* £ se 7¢ * = ee 


t 
f 
! 
: 


And 
Ra 2 Re (B4) 


where f is the propagation focal distance obtained in the usual manner from 
the SOQ code, and Z is the "new" scaled propagation distance. 


These changes are incorporated in the SOQ code primarily in subroutine 
QUAL, as documented by the following Fortran changes. 
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16. *ID ECSFIX 


The updates for ECSFIX are included to correct original errors in 
dimensioning Level II variables, and to reduce the resident array sizes at 
load and execution of the code. 
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*ID SEGSOQ 


The SOQ code, as currently configured, is too large to run on the Cyber i 
176 under AFWL Small Core Memory (SCM) restrictions (high speed core). The 
segmented load option of the CDC NOS/BE loader has been used to reduce : § 
execution time SCM requirements without loss of generality of the code. A 
segmentation loader, and the appropriate tree” structure of the code segmen- 


tation is required to take advantage of this feature. 


To incorporate this scheme into the SOQ code, the SOQ code required 
additional 3SLOBAL commons to‘save certain values, as described on the follow- 
ing Fortran listing. A segmentation tree was developed and is listed also. 
Further information on segmentation is available in the CDC/NOS/BE loader 
reference manual. This approach was selected instead of overlay structure 
tecause it is a more powerful tool, even though it is machine specific. 
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18. *WNDOW 


The aerodynamic window subroutine of GDL is used to model the effect of 
an aerodynamic window on the propagated field in a Monte-Carlo sense. The 
aerowindow subroutine simulates a random phase transmission function whose a 
rectangularly distributed random phase information can be selected with 
arbitrary "strength" or variance. This version of AEROW is designed to ae 
simulate the phase field degradation with rectangular probability distribu- 
tion in phase of 0.25). 
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